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Abstract 



This work is a continuation of [5]; it deals with rough boundaries in the simphfied 
' context of a Poisson equation. We impose Dirichlet boundary conditions on the periodic 

, microscopic perturbation of a flat edge on one side and natural homogeneous Neumann 

boundary conditions are applied on the inlet/outlet of the domain. To prevent oscillations 
on the Neumann-like boundaries, we introduce a microscopic vertical corrector defined 
in a rough quarter-plane. In [3] we studied a priori estimates in this setting; here we 
fully develop very weak estimates d la Necas [16] in the weighted Sobolev spaces on 
an unbounded domain. We obtain optimal estimates which improve those derived in 
[3]. We validate these results numerically, proving first order results for boundary layer 
^ ' approximation including the vertical correctors and a little less for the averaged wall-law 

, introduced in the literature [TH [T7] . 
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H ; 1 Introduction 

Cardio- vascular pathologies of the arterial wall represent a challenging area of investigation 
since they are one of the major cause of death in occidental countries. In this context, we are 
strongly interested in the accurate description of blood-flow characteristics in stented arteries. 
Specifically, we aim to understand the influence of a metallic wired stent (a medical device 
that cures some of these pathologies) on the circulatory system: our goal is to give a detailed 
description of the flow upward, inward and backward the region of stent's location. Actually 
the stent could be seen as a local perturbation of a smooth boundary of the flow field. The 
change, from perturbed to smooth, strongly contradicts the hypothesis of periodicity faced 
by the author in [H [5] . 

Although this problem was tackled in [T2] , our formalism follows ideas presented in [19] 
for interior homogenization problems, and it should be easy to extend it to other linear elliptic 
operators. A first step in this direction was made in [3] for a simplified Poisson problem: we 
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set up a formal approach to handle natural boundary conditions at the inlet and the outlet of 
a straight rough domain; then we proved rigorously, via specific a priori estimates, that the 
boundary layer approximation - built by adding some vertical correctors - converges to the 
exact solution of the rough problem. These estimates validated our approach. 

In im I12j . the authors introduced, via very weak solutions [TB], estimates of the 
error between various approximations and the exact rough solution. These estimates were 
established in a piecewise-smooth domain Q,^, limit of the rough geometry, when the roughness 
size e goes to zero. For a fixed e, this approach allows to estimate the error of an effective 
wall law approximation defined only in the smooth domain. In [3], we did not obtain optimal 
estimates in the norm. The major difficulty was some dual norm of a normal derivative 
as explained below. The present work fills this gap. 

In section [21 a short presentation of the problem and the material introduced in [3] are pre- 
sented. The difficulties that this paper overcomes are then faced in the next sections: firstly, 
the microscopic approximations live on unbounded domains and thus belong to weighted 
Sobolev spaces. As a result, one needs to derive very weak solutions on a quarter-plane, in 
these spaces (see section [3]) . Then one should connect these microscopic very weak estimates 
to the macroscopic problem we are really interested in. At this scale, the approximations live 
in the bounded domain $7^ and regular solutions belong to a specific subspace of H^{^^). 
While this correspondence was introduced for fractional test spaces in [3], here it is extended 
to the trace spaces specific to the regular solutions above (section |4]) . In section [U we analyse 
the convergence of the full boundary layer approximation towards the exact solution using 
arguments introduced in the previous sections; optimal estimates are obtained. Then for the 
first order wall-law, the convergence rate is shown to be equal to the one obtained in the pe- 
riodic case [3]. In a last part, we provide a numerical validation of the theoretical results. We 
compare various multi-scale approximations with a numerical solution of the complete rough 
problem. This comparison is made in Sobolev norms for various values of e. An accurate 
control of the mesh-size with respect to e and a P2 Lagrange finite element provide twofold 
results: the full boundary layer approximation shows the maximal convergence rate that one 
can expect from our numerical discretization, however, the standard averaged wall-law shows 
poorer results than expected. 

2 The framework 
2.1 The rough domain 

We set a straight horizontal domain 0,^, defined by 

:= G s.t. xi g]0 : 1[ and ef (^^^ < 3:2 < l} , (1) 

where / is a Lipschitz continuous function, 1-periodic. Moreover we suppose that / is bounded 
and negative definite, i.e. there exists a positive constant 5 such that 1 — 6 < f{yi) < S for 
all yi E [0, 27r]. The lateral boundaries are denoted by Fin and Fout, and their restrictions to 
]0, 1[, FinP (resp. F^^^). The rough bottom of the domain is called 

F^ :={xEM2 s.t. x2 = e/(^)}, 

while the top is smooth and denoted by F^. In the interior of the domain one sets the square 
piecewise-smooth domain Q,^ :=]0, whose lower interface is denoted by F*^, (see figH]). 
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Figure 1: (macroscopic) Rough, smooth and (microscopic) cell domains 



2.2 The exact rough problem 

In order to identify more precisely the influence of vertical non-periodic boundary conditions, 
we consider a singular perturbation of a linear profile: we look for solutions of the problem, 
find u G H^{Q.^) such that 

' - An' = 0, in 

< u^ = U, on r\n' = on T\ (2) 
^ d^u" = 0, on Tin U Tout- 

When e goes to zero we recover the linear profile v9 = 11x2, while this profile is explicit 
what follows shall apply with few modifications to the case of an implicit function which 
solves the problem: 

r - An° = 0, in 

< u° = [7, on = on ^^ 

. 5^u° = 0, on r;^^ ur^^f 

One can show (see [I2l Hj ) that 

Ih^ ~ ^°IL2(f70) ^ 

within the very weak solution framework d la Necas, that will be detailed below (see section 

2.3 First order approximation 

When one wants to improve the accuracy of the zero order approximation, one extends 
linearly using a Taylor formula in the neighbourhood of the fictitious interface T^. So we have 
= Ux2 for every x in Q*^. As the Dirichlet condition is no more satisfied on V", one should 
solve a microscopic problem that reads: find /?, whose Dirichlet norm is finite, such that 

- A/5 = 0, in Z+ U r U P, 

(5 = -2/2, on po, (3) 
/3 is yi — periodic . 

where Z+ :=]0, l[x]R+, P := {y g M? s.t. yi G]0, 1[, f{yi) < ?/2 < 0} and P° := {y G 
s.t. yi g]0, 1[, 7/2 = fiui)} (see fig. [1] right). In the literature this problem is widely 
studied (see [HI HZl Hj), so we only sum up the main properties of /?. 
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Lemma 2.1. There exists a unique solution of problem ([3]). Moreover, 



lim P{yi,y2) = P for every yi, and [3 := I P{yi,0)dyi. 

The convergence is exponential and one has a Fourier decomposition: 
oo „i 
m= E /3fce2-(-l^l^^+*^^^),VyGZ+, where (3u := / I3{yi,^)e^^''^y^ dy^. 
k=-oo -^0 

If u'' were periodic, we could set the first order approximation to be 



l,oo , 



(1) 



But this does not satisfy the homogeneous Neumann boundary conditions on Tin U Tout when 
approximating the solution of ([2]). In [3] we introduced a vertical corrector. We denote it by 
.^in; it solves the problem: 



- A^in = 0, in n, 
9i.6n(0,y2) = (0,2/2), on£^, 
, 6n = 0, on B. 



(5) 



where we set 11 := U^^[Z+ U F U P + kei]. The vertical boundary is denoted by E := {y ^ 
n,?/i = 0} and the bottom by B := U^^{y G P° ± kei} (cf. fig ED- In what follows we will 
write n' := M^, B' := M+ x {0} and E' := {0} x For the rest of the paper, we define the 
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Figure 2: Semi infinite microscopic domains: II, the rough quarter-plane and H' , the smooth 
one 



usual Sobolev space: 

W^'P{n) := {^v E V'{n) s.t. \D\\{1 + e LP{n), O < |A| < m } 

where p is a distance to the point (0,-1) exterior to the domain 11. Shifting the latter point to 
(0, 0) gives an equivalent norm so that we will not distinguish between these two distances. 
We refer to [9l [151 1] and references therein, for the detailed studies of the weighted Sobolev 
spaces in the context of elliptic operators. 

In the first part of this study [3], we have rigorously shown the results regarding ,^in: 
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Theorem 2.1. There exists a unique solution G Wc];'^(n) of problem ([5]) where a is such 



that \a\ < aQ := \/2/7r, 



moreover 



I6n(y)| < TTT^' e n s.t. p{y) > 1 

(l + p2(y)) 2(1-217) 

where M is a positive constant such that M < 1/(1 — 2a) ~ 10. 

This theorem is based on Poincare-Wirtinger estimates for the Sobolev part and a Green's 
representation formula in a quarter-plane 11'. Combining this two arguments, one obtains 
estimates on the decreasing properties of ^in. In the same way we define the vertical boundary 
layer corrector for Tout that solves the problem: 

-Aeout = 0, inn_, 

du(,out (0, y2) = - dyl3 (0, 2/2), on (6) 
,'?out = 0, onS_. 

where we set n„ = U;;^^[Z+ U T U P - /cei], the bottom being denoted by = U^^{y G 
— kei}. In Theorem, 12.11 as everywhere else in the rest of the paper, the properties 
derived for are equally valid for ^out- Thanks to these correctors, one completes the 
previous boundary layer approximation by writing: 



l + e(5J dx2 
This approximation satisfies the problem: 
An,^'°° = 0, on 0^ 



T,OUt 



Xi-l X2' 

e ' e 



(7) 
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In other words, we correct the 0(1) error on the normal derivatives of the boundary layer /3 
on Fin U Tout by the normal derivatives of the correctors at a distance yi = ^ of the vertical 
boundary E. On F^, the main errors are due to and ^outi the contribution of f5{x/e) — f5 
being exponentially small on F^. 

In [3] we set up adequate tools to handle a priori estimates for the error. We adapt them 
to the specific boundary conditions in problem ([2]), and claim 



Theorem 2.2. The boundary layer approximation ul''^ satisfies the error estimates in the 
Dirichlet norm: 



,l,00"i 



W^) < ke-'<'^-'l-^) 

where the constant k is independent of e and the constants a and M are defined as in Theorem 
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While in Theorem 5.1 [3], the convergence order is only e, here, it is improved because 
the perturbed profile is only linear: there are no second order errors in the sub-layer. The 
proof follows exactly the same ideas as in Theorem 5.1 in [3]. The estimates for the very weak 
solutions presented in Theorem 5.2 in [3] are not optimal and this is the main concern of the 
present article. 



3 Dirichlet-Poisson problem in a quarter-plane 

Domains, coordinates and notations We define the shifted domain 11^ ■=]l, cxo[xM+ and 
its boundary dU\ := E[ U B[ where E[ := {/} x M+ and B'l :=]/, oo[x{0}. We denote by Hj, 
the 7r/4 radians rotation of IIJ with respect to the origin (0,0) and 9nj, the corresponding 
boundary. M.f, := Mx]/',oo[ will represent the straightening of nj/, and I' := l/V^. The 
domain IT;, can be parametrised as ITj/ = {z ^M? s.t. 2i G M and Z2 > I' + \zi — l'\ =: a(zi)} 
whereas the change of variables from IT^, towards M^, is given by u; = W{z) s.t. 

j Wl = Zl, 

yW2 = Z2 — \zi — 

Later on we will also need the regularised version of domains above 



:= {yen; s.t. y2>^], 

z en;, s.t. Z2>1' + ^{zi - I'Y + s2 a,(zi)| 



(8) 



and the corresponding boundaries are set according to this definition. The mapping straight- 
ening nj, J, to Mj^, is set as w; = >Vs(z): 




3.1 Weak solutions 

We consider the solution of the problem: find v in Wo'^(n;) solving 

f A?; = 0, in rt, 
yv = g, onOT;, 

where (7 is a function belonging to wY^idX^i). We emphasise here that, as in Chap. 5 |16j . 

2 

we require a little more regularity on g than the usual fractional norm. As we work with 
weighted Sobolev spaces and we control the tangential derivatives of the data, the existence 
and uniqueness results for weak solutions are nor standard neither so straightforward. As 
they will be used extensively in what follows we provide a detailed presentation. 

In order to give a variational formulation of problem ([9|) , we need to construct lifts of the 
boundary data in the weighted Sobolev context. As we apply changes of variables above, we 
have to insure the compatibility of weights with respect to these mappings. Thus we present 
a detailed adaptation of the results for the half plane introduced in [9j . 

To solve problem ([9]) we use the Poincare-Wirtinger estimates because the weighted log- 
arithmic Hardy estimates are not valid in the specific case where a + n/p is an integer in 
W2'^(J1[) (see [HE] and references therein). 
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3.1.1 Weighted Sobolev extensions 

The first change of variables presented above, z{y), is a rotation of the domain around the 
origin: it preserves the distances to the origin. Once straightened in we are exactly in the 
position to construct extensions introduced by Hanouzet in Theorem II. 2 of [9], so we set: 

V{w)= [ K{t)g{t{w2-l') + wi)dt, 
J\t\<i 

where $ is a cut-off function such that 

Suppa>G [0,^[, $(0) = 1, $GC-([0,i]), 

and K is a regularising kernel i.e. K € C^(] — 1 : 1[) and K{s)ds = 1. Our lift then reads: 

R{g){w) =^i{w)V{w), V-wGM^i 

Then one has: 

Lemma 3.1. For every function g G VFi'^(R) one has: 

2 

\\R(9)\\w^^.^^ul) ^ ''Mwl'\mx{i'}V (10) 
where the constant k is independent on I' . 

Proof. According to the definition of the norm associated to W(}'^(nJ), 



/ (-) 



dw <k [ ( — ] [ g'^{t{w2 - I') + wi)dtdw 

mf, \ pj J\t\<i 

- iMx[o,i] 1 + ^? + (^0^ + + m ^ ' 

• / g'^{txJl + wl + {Uy + wi)dtdxdwi 
J\t\<i ^ 

< k — ^ dwi < k\\q\\,,rl,2,^a^, 



where we used the change of variables x = {w2 — I') / \/ 1 + wf + (/')^ and a shift wi := 
tx^yi + wf + {I'Y + wi in a second step as suggested in Lemma II.2 [9]. If 5 is a regular 
function, the gradient is estimated as: 



|V(W)| < 



\{V^)V\ 
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So the first integral proceeds as above, whereas in the second part, performing the same 
changes of variables, one gets: 



/ ■^'^\VV\dw<k [ [ g{t(:w2-l') + wi){l+tfdtdw 

JRl Jr^, J\t\<l 



\t\<l 

<k [ Jl + wl + {l'y<^^{x) [ (g' (txJ l + wl + {I'Y + \ dtdxdwi 
^Rx[0 ^ J\t\<l V V ^ // 

< k [ ^l + wl + {l'Y{g'fdwi < A;||5lUi,2(j,) 

JR. ^ 

extending this to VFi'^(M) functions by density arguments ends the proof. □ 

2 

In order to use these estimates in Yi[ we have to guarantee that they apply also when the 
domain is a quarter-plane. 

Lemma 3.2. For every g G VFi'^((9n^) there exists a lift denoted by R{g) in Wq''^(J1[) such 

2 

that 

where the constant k is independent on I. 

Proof. As mentioned above the rotation does not change the distances, so one has to consider 
the mapping from M^, to ITj, (resp. M x {I'} to dtl'i,) from both sides of (jlOp. The straightening 
by the continuous piecewise linear transform w = W{z) is defined in ([HI); one has that 

\V^R{g)\dw= [ {V,W^^V,W~'"V,R{g),V,R{g))det{V,W)dz. 

The eigenvalues of VzW~^VzW~"^ read A± = (3 ± ^/5)/2, they are positive definite and 
independent on Thus there exists a constant k such that 

On the other hand, one should focus on the equivalence of trace norms between wl''^{dlii,) 

2 

and WP(M x {/'}): on dU'i,, 1 + p^{z) = 1 + (zi)^ + (/' + \zi - This gives the existence 

2 

of a constant k independent on /' s.t. 

i + zj + {I'f < p{z) < k{i + zl + (I'f), vz G an;, 

In turn this implies that 

T + {l + p'{z))k9'Ada{z) 



en;, [(l + p2(2))2 



Jrx{1'} [{l+wj + {l')^)^ I 



□ 



Remark 3.1. The previous lemma applies with only minor changes to the case of the smooth 
domain sequence (IIJ 5)sg]o.i] • 
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3.1.2 A priori estimates 

At this point we are ready to prove the existence and uniqueness of a weak solution of problem 
Q. We denote by W^o'^(n9 the subspace of W^o'^lnj) such that: 

Wl^i^i) = {v^ W^o'^(n;) s.t. = on OT^}. 

Proposition 3.1. There exists a unique weak solution v G VFQ'^(n^), of the problem ([9]), 
moreover one has: 

where the constant k is independent of I. 

Proof. Using the lift of g given above, ([9]) becomes: find v G WQ^'^(n^) s.t. Av = —AR{g). 
Testing this equation hj tp £ I'(n^) one has that 

{Vv,V^) = -{VR{g),V^), 

by density of 25(11^) functions in W(j'^(nJ), the r.h.s. is a linear form on VFQ'^(n^) and the l.h.s. 
is a bi-linear bi-continuous form on the same functional space. Thanks to Poincare-Wirtinger 
estimates in a quarter-plane the semi- norm is actually equivalent to the VFQ'^(n;) norm. Thus 
one has the existence and uniqueness by the standard Lax-Milgram theorem. Moreover, one 
has that 

where all the constants do not depend on I. Subtracting R{g) one gets the desired result. □ 

3.1.3 Regularised problems 

In the rest of the paper, we need a little more regularity in order to construct a Green's 
formula adapted to the Lipschitz domain Ilj (see the very clear and detailed explanations of 
§ 1.5.3 in [S]). Thus in this paragraph, we construct regular approximations of problem Q. 
This is done by approximating IIJ by a sequence (11^ g)sg[o^i] of C°° domains defined above. 

For every given function g G Wl:\dn[,), one sets gs G Wi' (911; g) by writing 

2 2 ' 

gs{zi,as{zi)) := g{zi,a{zi)) = g{zi), Vzi G M 
and, where it is not ambiguous, we will drop the s and use g instead of gs- 

Lemma 3.3. The approximating sequence of data {gs)se[o,i] stable with respect to the 
W\''^{dl\'i,) norm: 

2 

whereas for the lifts one has 

Il^s(5s)llwi.2(n/ ) <k'\\g\\^i2^^Q^^y 

L .a ^ I 

and the sequence {Rs{gs))s converges to R{g) in the VFQ'^(nJ,) norm: 

Vry 35 > s.t. < s < 5 =^ \\Rs{gs) - R{g)\\y^i.2 (j^.^^-^ < V- 
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Proof. For s small enough, one has: 

a{zi) < as{zi) < 2a(zi), Vzi G M. 

Similarly it is easy to show that \a'g\ < \a'\. This gives the first claim. The proof of the second 
claim is identical the proof of Lemma [221 We extend Rs{gs) in HJ, \ ITJ, ^ by g. It is easy to 
prove that there exists a constant k s.t. 



\9\\w\-^{dfi')- 



Up to a subsequence, the compact imbedding of T4^Q^'^(nj,) into W'^i{Il'i,) implies strong 
convergence in the latter norm. The main focus is the semi-norm convergence. As seen in 
the lemma above the norms are consistent when passing from ITj, to M.f,, and the same holds 
when passing from 11^, ^ to M^, for the same reason. Making the change of variables Ws, one 
can re-express both lifts in R^, as: 

Rsig) = $ 



W2-1' 



+ ^1 + in 



f\2 



V{w), 



where Usiwi) = \/s'^ + wi — \ wi\. For the rest of the proof we set w = {wi,W2 +uJsiu;i)). We 
decompose I := Vw{Rs{g) — R{g)) in four pieces: 

4 

I :=Y^ Ij :=V^(^(w)) - '^{w))V{w) + V<i'{w){V{w) - V{w)) 

+ (^'(u}) - ^{w))VV{w) + '^{w)V{V{w) - V{w)). 
The first three terms can easily be estimated by s'''||g'||^i,2^jjj^ where 7 is some positive con- 

stant, thanks to techniques used in for fractional trace spaces. In our case those terms are 
even easier to treat because no fractional norm has to be used. The term /4 is more delicate 
because there is no derivative left, so we have to use the continuity of weighted translation 
operators; here we give the sketch of the proof: 



Ia<2 



"i/iw) / K{t){g'{tw2 + wi) - g'{tw2 + wi)} dt 
J\t\<i 



h.i + 1 



^(w) / K{t)g'{tw2 + wi)tJ^{wi)dt 
J\t\<i 

the term 14^2 is estimated again as the other terms, we focus on 14^1 

J := lldw<2k ^"^(x) \ g' [txJl + wl + iVf + uJs{wi)+wi\ 

Jul Jrx[o,^] I V ^ / 

-g' (^tx^l + wl + {l'y + wi^^ dt, 
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where, as before, we set x := {w2 — I') / \/ 1 + w\ + {I'Y ■ Then we use a version of Lemma 
II. 2. in [9] extended to all types of powers of the integration weight to conclude that: 

J <k I ^Jl + wl + {g{wi + tzuj{t, z, wi)) - g' {wi)Y dtdzdwi, 



where uj{t,z,wi) = uj{wi{t, z,wi)). At this point, we can follow the proof of Theorem 2.1.1 
in [16] that states the continuity of the translation operator in L^. The only difference is that 
the term tzoj{t, z,wi)) itself depends on the integration variables. This dependence problem 
is overcame by noting that 



Vr/ > 0, 35 > s.t. \tzCj{t, z,wi))\ < ^< 6, V(t, z, wi) e] - 1, l[x 
this in turn implies that on the set of continuity points of g' , one has 



1 



\g'{wi + tzuj{t,z,wi)) - g'{wi)\ < -, 
and the rest follows exactly as in Theorem 2.1.1 in [16j . □ 



3.1.4 Convergence 



Now we state the existence and uniqueness of the regularised problem: find Vg G VFo^'^(5nj ^) 
s.t. 

f Avs = 0, in n; 

\ (11) 
[vs= g, on dHi ^. 

Proposition 3.2. For every fixed s G [0, 1] there exists a unique weak solution Vs G VFQ^'^(n^ ^) 
of problem (llip . satisfying 

V?7 > 36 > s.t. s < 6 =^ \\vs — f ||^i,2^j-j,^ < r], 

where Vg is extended by g in 11^ \ 11^ ^ . 

The proof is a straightforward consequence of Proposition 13.11 applied to the regular do- 
main g and of Lemma 13.31 for the convergence part. 

If we now restrict ourselves to the case where g G £(Ili), this latter space being dense in 
VFi'^(nJ) (see for instance Theorem I.l in [9]), we get more regularity, namely: 



Lemma 3.1. If g G S{dll'i) then Vg, the unique solution of problem (jlip . belongs to H^^^ili'^ ^ 
for every fixed s g]0, 1]. 

Proof. It is easy to show that if 5 G £{dYl'i)^ then Rs{g) G £(J^i) and thus by standard interior 
regularity results one gets that Vg ■= Vg — Rsid), which belongs to li^Q^'^(nj j,), is actually in 
^loci^i s) [3 El US]. Moreover because of the Dirichlet condition and the C°° regularity of the 
boundary, the regularity of Vg can be extended up to the boundary by the same method. □ 
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3.2 Weighted Rellich estimates 

Above, we constructed the tools necessary to adapt the very weak solutions presented in 
Chap. 5 in [16], to the weighted context. 



Proposition 3.3. Let s G]0, 1] be fixed. If g G '^(11^^) then Vs G -ffjQ^(nj^) and one h 



as 



moreover 



II 9uVs Iliy0.2(9n; J - ^lls'llvKi'^(9nj j- 

The operator T defined from £{Il[ J as T{g) = dyVg is extended by continuity into a mapping 
from Wi'^(5n; J on Wl''^{dU[J and one has that 

2 ' 2 ' 

ll^(5)llH/?-2(an; J ^ f^Mwl^^dui 
where the constants k do not depend on I. 

Proof. We rotate again IIJ ^ by 7r/4 radians to switch to the chart (zi,Z2)- the boundary of 
dll'i g is expressed as {zi,as{zi)). We set the partition of unity, 

+ 00 

^iPr{z) = i, yz€ui, 

r=0 

the functions (fr being defined as: 

:= e yp-^'-^J l[2r-i,2'-](p) + l[2^2'■+l](p) + e ^'^ ^ l[2r+\2^+2]{p), r > 1, 
V'o := l[o,2](p) + e VP-*; l[2,4](p), (12) 
fr ■■= ——, , Vr > 1, (fo: 



j=r-l 

where by 1 5 we denote the characteristic function of a given set S. Then we define 

+00 



h^ := (0, —ipr{z)p{z)), and h := /i''. 



r=0 



Thanks to Proposition [XH Vs G ^^^^(n^, J; one is allowed to set locally the Rellich formula 
([16j. p. 245). When adapted to the Laplace operator, it reads: 

((/i*" ■ i^)ld2 - (h® + 1^ <^ h)Vvs,Vvs) da{z) 
= [ (div h"" Ida - (V/i'' + {Vh''f)Vvs,Vvs) dz + / (/i^ • Vt;,) At;, dz. 
We have that 

{v, h') = ifr{z){l + (a'j2)-|p(^) > A:(^^p(z) (13) 
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because \a'g\ < 1. Developing the boundary term in normal (z^) and tangent (r) directions, 
one has that, in fact. 



{K" ■ v){d,vsf + 2{K' ■ T)drVsd,Vs - [h^- ■ iy){d,vsf da{z) 

s 

= - [ (div Ids - (V/i'' + {Vh''f)Vvs,Vvs) dz. 

The first term in the r.h.s. above is estimated from below thanks to (|13p . the second and the 
third ones by their absolute value, giving: 



ipr{z)p{z){duVs) da{z) < I 2iprp{z){\drVs\\dyVs\ + {d^Vs) ) d(T{z) 

- [ {{div h''ld2 - (Vh' + {Vh''f)}Vvs,Vvs)dz. 

Then we sum with respect to r and apply the Beppo-Levi theorem for the boundary terms. 
For the interior r.h.s. above, due to the specific choice of cut-of function in ()12|) . one can pass 
to the limit with respect to the summation index r applying the Lebesgues theorem. These 
justifications allow us to write 

/ p{z){d^Vsfda{z)<f 2p{z){\d^Vs\\drVs\ + {d^Vsf) da{z) 



L ({d-(")M.-(V + V-)(0)}v..,V..).., 



note that it is important here that the derivatives of h contain only p but no cut-of function, 
this explains why we don't estimate the interior terms before summing over r. By Cauchy- 
Schwartz one gets that 

II duVs \\w°-^{dn[ J - ^ I II^'IIh/i'^cot; jH ^"'"^ \\w°'^(du[ j + ll^«llH/o'^(n;, j 

Thanks to the Young inequality and Proposition 13. 1|. one obtains the first estimate of the 
claim. Because the r.h.s. only depends on the wl''^{dll'i g)-norm, the result can be extended 

to every function g belonging to Wi' {dH'i ^) by density and continuity. □ 

Remark 3.2. Note that this result holds in fact for any polynomial of p and what follows 
could be extended as well to any weighted Sobolev space. One only needs to choose the proper 
scaling for the cut-off functions ipr with respect to the weight. 

Proposition 3.4. // (vs)se[o,i] sequence of solutions of problems (jlip . then the next 

properties hold: 

(i) There exists a constant k dependent neither on s nor on I such that 

II duVs \\w°''^{au[ J - ^II5'IIw'J'^(ot;)' 
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(ii) there exists w G Wi (SIIJ^) s.t. 

2 ' 

d^Vs ■ccr, in W^'^{dl\'i g), 

2 ' 

(Hi) for every function u G Wo^'^(nJ) one has at the limit s = 0, the Green's formula: 

/ Vv-'Vudy= / zuuda{y), 
Jn[ Jdu[ 

where v is the solution of problem ([9]). 

Proof. Part (i) comes from Proposition [3i3] combined with Lemma [3.3i Again we approximate 
g G VFi'^((9n^) by G £{dll'i) and we call the unique solution of problem (jlip with data 

2 

g^ given on By continuity of the solution of problem pT]) with respect to the data, one 
easily shows that 



V{vl-vs) 



< k 



Rs{g') - Rs{g) 



< k' 



1 



< k" 



g'-g 



W\'^{dll[) 



thanks to the weighted Rellich estimates above, one has also that 



7 



< k 



g'-g 



And because g^ is regular, G i?jQj,(nJ^). This allows us to write the Green's formula for 
every h G V(J\\^^): 

[ VvlVhdy= [ dyvihda{y) 
Thanks to strong convergence shown above, one can let 5 go to the limit and get: 



VvgVh dy 



en; 



duVshda{y) 



(14) 



note that working only with (511^ one can not write directly the Green's formula: 

2 ' 

is not regular enough. Thanks to the last estimate of Lemma 13.31 one has that 



lim / Vvs^hdy 



VvVhdy 



Thus a limit for the boundary term in the r.h.s. of (I14p exists. We express the boundary 
term for a fixed s in ^ in z coordinates, and we choose h to be only a function of zi on the 
boundary. Moreover u := (sgn(2:i), —1) is the limit outward normal, (resp. Vg ■= {a'(zi), — 1)) 
and we write 



lim / Vv{zi,as{zi)) ■ i>h{zi) dzi = / VvVhdy + 



+ lim / \/v{zi,as{zi)) ■ {D - Ds}h{zi) dzi. 
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The last term can be estimated through the Cauchy-Schwartz inequahty 

Vv{zi,asizi)) ■ {D - Ds}h{zi) dzi < ||5f||^i,2(gn;) || (sgn - "■s)Hw°iidU[) 

Applying the Lebesgues theorem, the last term in the latter r.h.s. goes to zero. Indeed, note 
that, by Poincare-Wirtinger arguments, there exists a constant k s.t. 



\h 



By density and continuity arguments, one extends the Green's formula to all functions in 
W^'\U[). □ 

Proposition 3.5. Let T be the linear continuous operator from wl''^{dlli) on Wi''^{dll'i) 

2 2 

S.t. T{g) = duV where the normal derivative is to he understood as a weak limit exhibited 

2 12 12 

above. Then T is extended as a map from W \ {dli[) on W~i' (9nj) where W~i (SIIJ) = 

{w\^\^^)y. 

2 

\ '2i 2 12 

Proof. Wi' {dli[) is dense in W^^{dli'i). Let g,h be two functions of W"i' (511^) such that 

2 2 2 

f An = 0, in H;, f = 0, in Hj, 

< , and < , 

\u = g, on aHj, \v = h, on all;. 

Then the Green's formula from Proposition 13.41 applies twice, giving 

u{dyv)da{y) = / {d„u)v da{y). 

JdU[ 

Thanks to the Rellich estimates, one then gets 

{dyu)hda{y) < k\\9\\w°-l{du'.)\\^\\wl'\dn',)^ 
which gives that 

\\Tig)\\w-l'\du[) ^ ^\\9\\w°'Udn[r 
and density arguments complete the proof. □ 

3.3 Weighted dual estimates on the normal derivatives 

We now return to the study of the vertical boundary layer corrector that solves problem ([5]). 
Thanks to Proposition 13.51 we derive one of the key point estimates of the paper: 

Proposition 3.6. There exists a unique solution ,^in G TyQ'^(n), of problem ([5]). Moreover 
there exists a constant k that does not depend on I s.t. 

-I i_ 

II du^in lliy-p(OT;) - ^ ( 7 



where the constant M is defined as in Theorem \2.1 
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Proof. The proof is a straightforward apphcation of Proposition 13.51 

II f^t'Cin ||,4/-i,2(gn,) < ^il6n|lty'''2(OTJ)- 

As is at least inside the domain, we use the point-wise L°° estimates from Theorem 
12.11 which give: 

r+oo ^ 

That provides the desired result. □ 

Remark 3.3. This result express the decrease of the normal derivative of on a vertical 
interface located at yi = I. These estimates improve the convergence rate obtained in Propo- 
sition 4 in J3i by a factor of almost y/l/l. Indeed we consider here the W l''^{dll'i) norm, 

~ 2 

--,2 

while in J^, only the Wq ^' (dH'i) norm was used. In the rest of the article we imbed and 
exploit the result above into the macroscopic very weak setting. 

4 Correspondence between macro and micro Sobolev norms 

In [3], a correspondence was shown between Hq (Tin U Tout) and a subspace of W^' {dlli), we 
extend it here between HQ{r[^ U T'^^^) and a subspace of VFi'^(5nJ) test functions. In what 

2 

follows the same could be written for T'^^^. Taking v G i?g(rin') we set 
V = V ("'2^2^ = V{0,X2), VX2 G [0, 1], 

and we extend v by zero on 911^. Note that this makes sense because v is zero at X2 = and 
X2 = 1, so that one has 

Lemma 4.1. For a given function v £ i?Q (Tin') and v defined above, the following equivalence 
of the Sobolev trace norms occurs: 

Mwl'\dn[) ^ ^ll^llHi(rin') ^ f^'Mwl'^an'^) 

where the constants k,k' do not depend on e. 

Proof. We start from the macroscopic side, the other way follows the same. 

1 



/ v'^{0,X2)dx2 = e [ v^{0,ey2)dy2 = e [ ( - , 2/2 ) c^y: 
JVin Jo Jo \^ / 



^7 2 



where the constant k is obviously independent on e. Owing that dy^v = edx^v, the derivative 
part is shown similarly. □ 
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5 Very weak estimates for boundary layer and wall law ap- 
proximations 

Turning again to the macroscopic error estimates, one defines the error rl'°° := u'' — ul'°° 
where is the exact solution of problem ([2]) and ul'°° the boundary layer approximation 
proposed in ([7]). It satisfies the set of equations: 

' Arl'°^ = 0, in 
rl'^ = 0, on 

= -||(-^"»' (^,^) (^.;)) . on r'. (15) 

In order to improve L^(0*^) estimates obtained in [3], we use the material above to prove the 
main result of this paper: 

Theorem 5.1. There exists a unique solution re'°° € H^{^'') of problem (jlSp ." it satisfies the 
estimate: 



^^mm(|+a,2-5lj-) 



the constants a and M being defined in Theorem \2.1\ 

Proof. For any given function F G Li^iyt^), we solve the regular problem: find v S H\){^^) := 
{u G H^{Vl^) s.t. u = on r° U F^} such that 

- At; = F, in 17°, 

= 0, on r-^uro^t, 
, u = 0, on r°ur^ 

According to Theorem 4.3.1.4, p. 198 [S], u G H'^{n^)r\Hl,{Vl^) so that v G H^{dVP) and, 
thanks to boundary conditions on F*' U F"'^, v G //q (Fin U Fout)- We are now in the position to 
apply the Chapter 5 of [TB] to write that: 



j^^ rl'^Fdx = - (r^-, 5,t;)^ouri + ( ^> 



where by the brackets we denote the duality pairing H~^, i/g (FjnUFout) and by the parentheses 
we denote the scalar product in L^(F''uF^). By standard interior regularity results one easily 
gets that G H^^^{IV) (resp. ^out G ^ioc(-'^)) normal derivatives 

dudn GL2(0,1), ("resp. 9,Ut & L\0,1) 
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for every fixed e. Thus the duahty pairing becomes an integral: 

(d,r}'^,v) 



[ duiin \-, — ]v{x)da{x) - [ dui, 



out ( ~' ~ ) v{x)da{x) 







1 



-,y2 Vin -,y2 



1 



1 



-,y2 Vont -,y2 



1 



dy2 



< e 

< e 



~2 7 



^^OUt||,^1.2(gn'j) + II 



du£,'m \\w~'l'^{dn\) ~^ II ^'^'^out 



out lliy-i.2(an'j)l 

~^ 7 



\V\\ TTitr' 



2 7 



where and -Uout are the microscopic test functions associated to the trace ofv on TIj^UPq^^ as 
in section m One then concludes this part setting ? = 1/e in Proposition 13.61 The L^(r''ur^) 
scalar product has been estimated in [3] , using a priori estimates for the part whereas the 
part uses again L°° estimates from Theorem 12.11 □ 

A direct consequence of this result is 

Theorem 5.2. The first order wall law solving 

' Au^ = 0, in 0°, 



on r 



ef^, on r°, 

OX2 



(16) 



dyU^ = 0, on Tin U Tout, 



satisfies the error estimate 



.11 



- 



where the constant k is independent on e. 

The proof follows exactly the same line as in Theorem 5.3 in [3], but the result is improved 
thanks to the Theorem 15.11 above . 



6 Numerical evidence 

We define the rough bottom of the domain by setting / in ([1]) as: 

/(yi) = -l + ^sin(27ryi), VyiG[0,l]. 

This is obviously a Lipschitz smooth function compatible with the hypotheses of the claims. 
In what follows we look for a numerical validation of theoretical convergence results above: 
we compute for every fixed e G [0, 1] 

- uj^ a numerical approximation of u"^ solving a discrete counterpart of problem ([2]). 

- u^'^i^, the periodic full boundary layer approximation (it does not contain any vertical 
corrector) defined in (jl]) 
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- u^''^ , the full boundary layer approximation including vertical correctors defined in ([7]) 

- uj^, the averaged wall-law presented in (jl6p . and the zero order approximation. 

We use the finite element method code freefem++ [10], in order to compute n^, 6n,/i 
and ^out,/i- The P2 Lagrange finite elements interpolation is chosen. 

Microscopic correctors As f3, and ,^out are defined on infinite domains, we have to 
truncate these and set up proper boundary conditions on the corresponding new boundaries. 
For (3, this was analysed in [13] so that we only need to solve 

' - ApL = 0, in Z+ U r U P n {y G M^y2 < L}, 
(5l = -?/2, on P\ 

< 

Pl is yi — periodic , 
_ OuPl = 0, on {y2 = L}. 

The approximation Pl is exponentially close to P with respect to L in the Dirichlet norm (see 
Proposition 4.2 [I^)- For the vertical correctors we set the domain 11^ := Iln [— 1,L]^ (resp. 
:= n_ n [—L, l]^ and we solve the problem 

' - ACt = 0, in n, 

5.Cin(0,y2) = -5,/3 (0,2/2), onE, 
' Cin = 0, onB, 

^ d^Ctn = 0, on {yi = L} U {7/2 = L} =: Gl, 

the symmetric problem for being omitted. By Proposition 4 in |3] and Proposition 13.61 
above, one easily deduces the convergence result: 

Proposition 6.1. There exists a unique solution ^(^ G WQ^'^(n^) solution of problem (fT7|) . 
moreover one has 

W^tn - ?m||H/i'2(ni) < kL""", 11^4 - ^inllvi/Of (n^') - kL~'^'^2M 

where the constants k,k' are independent of L and a and M are defined as in Theorem \2.1{ 
' is the restriction of Yi^ to M+ x M+ . 

In figures [6] and [HI we display the meshes obtained after adaptative procedure, described 
below, for P^ and S,^^^. The total number of vertices used in the meshes for discretising 

and Cout,h ^"^^ 39000, 78000 and 79000. In the simulation of the horizontal top is 
set to L := 10. For ^ and the vertical interface is set to L := 20. The contours of 

corresponding solutions Pj^ and ^^^^^ ^ are displayed in figures [3 [9] and [lOl whereas the normal 
derivative dyP^^ and —d,y£,^i^ are shown to coincide along {0} x [— 1,1] in figure [3l 

We perform a single microscopic computation. Then we re-scale the boundary layer to 
the macroscopic domain setting 

Pthix) := Pl: (^) , Ctn,e,h{^) := (f ) , CL,e,k{x) := CL,h (f ) , Vx G 
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0.2 0.4 0.6 O.f 



Fi gure 3: Normal derivatives di,f5^ and dy^^^^ on the vertical interface 



We quantify the interpolation error with respect to e. 



L2(no) 



where s is a constant dependent on the boundary's regularity, and hm a fixed maximum mesh 
size on the microscopic level, independent on e. In the same way one can set 



out ) 



L2(n0) 



— ^ (^Ihoutj/i Cout ||j^2(nL') + ll^out ^out|lL2(n-^) 

— ^£ (^m||'?out||j:f2,i.(nl-) + ^-f'llCout ~ 



out||^^o,2j.j^i,^ ) < ^e, 



where u \s a, real parameter depending on the angle of the corner of n_ at (0,/(0)), and 
H^''^ the weighted space defined p. 388 Definition 8.4.1.1 [8], that takes into account the 
corner singularity of second derivatives of These estimates give an upper bound on the 



convergence rate for the full boundary layer Ue'°°, namely: 



l,oo 



l,oo I 



+ 



where is a macroscopic mesh size presented in the next paragraph 



(18) 



Rough solutions When computing numerical approximations of u'^, one has to play with 
3 concepts that are interdependent: h the mesh-size, e the roughness size, and corner singu- 
larities that depend on the shape of the domain. 

In the periodic case considered in [5], and for / G C°°(]0, 1]), in order to avoid that the 
roughness size goes under the mesh-size, one could discretise the solution on a mesh such 
that h < ce. Due to estimates on the interpolation error and ff^(O^) regularity, one obtains 
a good numerical agreement for convergence rates between theoretical and numerical results 
(see m- 

In the non-periodic setting, corner singularities occur near Tout • I^ order to obtain con- 
vergent numerical approximations of near Tout, one should refine the mesh in the neigh- 
bourhood of (1, e/(l/e)). At the same time, in the regular zones, the mesh-size should stil be 
refined at least linearly with respect to e (as in the peridic setting [^). This complicates the 
local size of elements with respect to the size of the mesh ([8] p. 384). Thus, simply setting 
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uniformly h := ce does not provide accurate convergence results. On the other hand, one aims 
to have a strong control on the mesh size far from the corner: for instance in these zones, the 
mesh-size could be fixed on a uniform grid. These considerations led us to use an overlapping 
Schwartz algorithm [18]; we split in two parts: ^l^ is discretised with a structured grid of 
size H := ke'^ (7 is discussed later), whereas a second domain reads 



r^' n <^ X G M s.t. X2 < 



and contains the rough sub-layer. On O^'*^ we perform mesh adaptation in order to capture 
geometrical and corner singularities. The maximum/minimum mesh-sizes are set: 

hrain = niin hx, /imax = min hK 

where hx is the diameter of triangle K in the triangulation Tk of Q^'^. At each step m of 
the Schwartz algorithm, we solve two problems. We set to be the solution of 

( = 0, in 17°, 

= 1, onr\ 

^y(m-l)^ On T", 



0, on r;„ u r' 



out J 



and Vf* solves 



AVr = 0, in 
= 0, on r% 

U^, on {x2 



VT=UT, on {^2 = 3^}, 
d,VT =0, on (Fin U Tout) n{x G 



s.t. X2 < 



To) 



and we iterate the procedure until 



{U^ - V^yda{x) < tol. 



(0,l)x{0}U(0,l)x{^} 

where tol is a constant set to 10~^°. During this step both meshes are kept fixed. 

Then we refine the sub-layer mesh Tk in order to account the corner singularity. This step 
provides a new mesh-size distribution updating /imin and /imax- We use adaptative techniques 
presented p. 92 of the freefem++ reference manual [10]. This procedure is compatible with 
the mesh requirements displayed in Theorem 8.4.1.6 p. 392 in [8] and guarantees standard 
interpolation errors with respect to the mesh size. 

We iterate these two steps: solve the Schwartz domain decomposition problem and then 
adapt the mesh. The iterative algorithm stops when /imax < H. Through this algorithm we 
insure both a given mesh size H and a refined mesh near the corner. 

We tested different values of 7 where setting H = ke^ , /c a is given constant, choosing 
7 > I does no more change convergence results below. We plot in fig. HI /imax and /imin as 
functions of e. The adaptative process gives approximately /imin ~ ce 



2.29 
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Figure 4: Mesh sizes /imin and /imax as functions of e 




Figure 5: Errors in the L'^{Q^h) (left) and the (right) norms with respect to e 

We plot in fig. [TTl the meshes obtained thanks to our iterative scheme for e G {^jI}- 
In fig. [121 we display the corresponding solutions u^. Next, we construct boundary layers 
using microscopic correctors above. We compute the errors u\ — , u\ — v} , u^^ — u^'^j^ and 

uf^ — u^''^ in the L'^{Q^) norms, and display them as a function of e in figlSl The numerical 
convergence rate, obtained by interpolating results above as a powers of e, is displayed in 
table m 



norm / approx. 




< - 


% - K*,h 






0.78783 


1.11 


1.1 


1.462 




0.787 


0.6869 


0.70 


1.346347 



Table 1: The errors convergence rates displayed as powers of e 



Discussion When the vertical correctors are not present, the boundary layer approximation 
is not only less accurate but also the rate of convergence is less than first order, the difference 
is visible in L'^{Q^) but is significant in the H^{Q^) norm. Nevertheless, and as explained 
above, when using a single microscopic computation of the correctors for every e, it is not 
possible to get better convergence results than ea. This is actually what we obtain for our 
more accurate approximation n^')^. This validates our theoretical results. The surprising 
phenomenon that we are at this point not able to justify is the poor convergence rate of the 
wall law , that should according to our estimates be e2. Observed in [5], performs even 
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worse convergence rate than in the H^{Q^) norm. The results of Theorem 12.21 are fairly 
approximated for what concerns the H^{Q^) error of u^''^ ■ 

7 Conclusion 

Our approach provides an almost complete understanding of the non-periodic case for lateral 
homogeneous Neumann boundary conditions in the straight case, (no curvature effects of the 
rough boundary |17j). A forthcoming paper should adapt these results to the case mentioned 
in the introduction: a smooth boundary forward and backward the rough domain via domain 
decomposition techniques. Another extension to the Stokes system should follow as well. 



Acknowledgements. 

The author would like to thank C. Amrouche for his advises and support, S. Nazarov for 
fruitful discussions and clarifications, and H. Teismann for proof reading. This research was 
partially funded by Cardiati^, an industrial partner designing and commercializing metallic 
wired stents. 



References 

[1] C. Amrouche, V. Girault, and J. Giroire. Weighted sobolev spaces and laplace's equation 
in M". Journal des Mathematiques Pures et Appliquees, 73:579-606, January 1994. 

[2] C. Amrouche, V. Girault, and J. Giroire. Dirichlet and neumann exterior problems for 
the n-dimensional laplace operator an approach in weighted sobolev spaces. Journal des 
Mathematiques Pures et Appliquees, 76:55-81(27), January 1997. 

[3] E. Bonnetier, D. Bresch, and V. Milisic. Blood flow modelling in stented arteries: new 
convergence results of first order boundary layers and wall-laws for a rough neumann- 
laplace problem, submitted. 

[4] D. Bresch and V. Milisic. High order multi-scale wall laws : part i, the periodic case, 
accepted for publication in Quart. Appl. Math. 2008. 

[5] D. Bresch and V. Milisic. Towards implicit multi-scale wall laws, accepted for publication 
in C. R. Acad. Sciences, Serie Mathematiques, 2008. 

[6] L. C. Evans. Partial differential equations, volume 19 of Graduate Studies in Mathemat- 
ics. American Mathematical Society, Providence, RI, 1998. 

[7] D. Gilbarg and N. S. Trudinger. Elliptic partial differential equations of second order. 
Classics in Mathematics. Springer- Verlag, Berlin, 2001. Reprint of the 1998 edition. 

[8] P. Grisvard. Elliptic problems in nonsmooth domains, volume 24 of Monographs and 
Studies in Mathematics. Pitman (Advanced Publishing Program), Boston, MA, 1985. 

[9] B. Hanouzet. Espaces de Sobolev avec poids application au probleme de Dirichlet dans 
un demi espace. Rend. Sem. Mat. Univ. Padova, 46:227-272, 1971. 

^ www. cardiatis.com 



23 



[10] F. Hecht, O. Pironneau, A. Le Hyaric, and Ohtsuka K. Freefem++. Laboratoire Jacques- 
Louis Lions, Universite Pierre et Marie Curie, Paris, 2005. 

[11] W. Jager and A. Mikelic. On the interface boundary condition of Beavers, Joseph, and 
Saffman. SIAM J. Appl. Math., 60(4):1111-1127, 2000. 

[12] W. Jager and A. Mikehc. On the roughness-induced effective boundary condition for an 
incompressible viscous flow. J. DiJJ. Equa., 170:96-122, 2001. 

[13] W. Jager, A. Mikelic, and N. Neuss. Asymptotic analysis of the laminar viscous flow 
over a porous bed. SIAM J. Sci. Comput, 22(6):2006-2028, 2001. 

[14] Willi Jager and Andro Mikelic. On the boundary conditions at the contact interface 
between a porous medium and a free fluid. Ann. Scuola Norm. Sup. Pisa CI. Sci. (4), 
23(3):403-465, 1996. 

[15] A. Kufner. Weighted Soholev spaces. BSB B. G. Teubner Verlagsgesellschaft, teubner- 
texte zur mathematik edition, 1980. 

[16] J. Necas. Les methodes directes en theorie des equations elliptiques. Masson et Cie, 
Editeurs, Paris, 1967. 

[17] N. Neuss, M. Neuss-Radu, and A. Mikelic. Effective laws for the poisson equation on 
domains with curved oscillating boundaries. Applicable Analysis, 85:479-502, 2006. 

[18] A. Quarteroni and A. Valli. Domain decomposition methods for partial differential equa- 
tions. Numerical Mathematics and Scientific Computation. Oxford Science Publication, 
Oxford, 1999. 

[19] E. Sanchez-Palencia and A. Zaoui. Homogenization techniques for composite media, 
volume 272 of Lecture Notes in Physics. Springer- Ver lag. , 1987. 



24 



Figure 6: The microscopic periodic cell after adaptative mesh refinement 




Figure 7: The microscopic periodic cell corrector (3^ 
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Figure 8: The microscopic domain of ^^yj after adaptative mesh refinement 



10.0168863 
■0.0236408 
10.0303953 
■0.0371498 
■0.0439043 
■0.0506588 
10.0574133 
■0.0641678 

■ 0.0709224 
10.0776769 

■ 0.0844314 

■ 0.0911859 
10.0979404 

■ 0.104695 

■ 0.111449 
10.118204 

■ 0.124958 

■ 0.131713 




Figure 9: The microscopic corrector .^^it 
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Figure 12: The rough solution computed for e G {g, 1/3} using a decomposition method 
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